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Abstract 

A new method is explored to detect extensive air showers: the measurement of radio waves emitted during the propagation of the 
5^ electromagnetic shower component in the magnetic field of the Earth. Recent results of the pioneering experiment LOPES are 
discussed. It registers radio signals in the frequency range between 40 and 80 MHz. The intensity of the measured radio emission 
is investigated as a function of different shower parameters, such as shower energy, angle of incidence, and distance to shower axis. 
In addition, new antenna types are developed in the framework of LOPES STAR and new methods are explored to realize a radio 
self-trigger algorithm in real time. 

Key words: cosmic rays, air shower, radio emission, radio detection 
PACS: 96.50.S-, 96.50.sd 



1. Introduction their acceleration and propagation mechanisms [1 2 3j. At 

energies exceeding 10 14 eV cosmic rays are usually studied 
An intense branch of astroparticle physics is the study h Y indirect measurements — the investigation of extensive 

of high-energy cosmic rays to reveal their origin, as well as air showers initiated by cosmic particles in the atmosphere. 

Different techniques are applied, like the measurements of 
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Fig. 1. Measured field strength as a function of the distance to shower 
axis for an individual shower |4]. 



particle densities and energies at ground level, or the obser- 
vation of Cerenkov and fluorescence light. An alternative 
technique has been recently revitalized — the detection of 
radio emission from extensive air showers at energies ex- 
ceeding 10 16 eV. 

Radio emission from air showers was experimentally dis- 
covered in 1965 at a frequency of 44 MHz [5J. The early 
activities in the 1960s and 1970s are summarized in [6]. 
Only recently, fast analog-to-digital converters and mod- 
ern computer technology made a clear detection of ra- 
dio emission from air showers possible. LOPES, a LO- 
FAR Prototype Station had shown that radio emission 
from air showers can be detected even in an environment 
with relatively strong radio frequency interference (RFI) 
[7]. Further investigations of the radio emission followed 
with LOPES |8|9ll0JllJ12j and the CODALEMA experi- 
ment |13|14| . paving the way for this new detection tech- 
nique, see also [15]. 

The LOPES experiment registers radio signals in the fre- 
quency range from 40 to 80 MHz [16J. In this band are 
few strong man made radio transmitters only, the emis- 
sion from air showers is still strong (it decreases with fre- 
quency), and background emission from the Galactic plane 
is still low. An active short dipole has been chosen as an- 
tenna. An inverted V-shaped dipole is positioned about 
1 /4 of the shortest wavelength above an aluminum ground 
plate. In this way a broad directional beam pattern is ob- 
tained. LOPES comprises 30 antennas [17] located on site 
of the KASCADE-Grande experiment |18|19j one of the 
best air shower experiments operating in the energy range 
between 10 14 and 10 18 eV. The LOPES data acquisition is 
triggered by large air showers registered with KASCADE- 
Grande. The latter measures the showers simultaneously 
to LOPES and delivers precise information on the shower 
parameters, such as shower energy as well as position and 
inclination of the shower axis. All antennas, including the 
complete analog electronics chain, have been individually 
calibrated with a reference radio source [20] . 

Most likely, the dominant emission mechanism of the 
radio waves in the atmosphere is radiation due to the de- 



Fig. 2. Scale parameter Rn as a function of the zenith angle times 
the mean distance to the shower axis 0]. 

flection of charged particles in the Earth's magnetic field 
(geosynchrotron radiation) |21|22|23] . In the frequency 
range of interest the wavelength of the radiation is large 
compared to the size of the emission region: the typical 
thickness of the air shower disc is about 1 to 2 m only. 
Thus, coherent emission is expected which yields relatively 
strong signals at ground level. 

2. Radio signal and shower parameters 

One of the primary objectives of LOPES is to investi- 
gate the measured radio signal as a function of parame- 
ters characterizing the extensive air shower. For this pur- 
pose the shower parameters measured simultaneously with 
KASCADE-Grande are irreplaceable. An empirical rela- 
tion has been found to express the expected east- west com- 
ponent of the field strength at a distance R from the shower 
axis as a function of shower parameters [24] 
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a is the angle between the shower axis and the direction of 
the Earth magnetic field (geomagnetic angle), the zenith 
angle of the shower, and En the energy of the shower induc- 
ing primary particle. It is interesting to note the absolute 
values of some parameters: the exponential fall-off has a 
characteristic length of about 240 m, much larger than the 
classical Moliere radius of electrons in air (^ 80 m). This 
indicates that the lateral distribution of the radio compo- 
nent is flatter as compared to the electromagnetic compo- 
nent, an important fact to build large-scale radio arrays 
with an economic antenna density. The measured radio sig- 
nal is almost directly proportional to the shower energy 
([0.95 ± 0.04] « 1). Such a behavior is expected for a co- 
herent emission of the radio waves from the air showers. 
This calibration of the measured radio signal is a first im- 
portant step towards the application of the radio detection 
as independent method to register extensive air showers. 
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Fig. 3. Sky map of air showers detected in east- west (left) and north-south (right) polarization direction [25]. The boxes mark the arrival 
direction of the air showers with the strongest radio signal. 



Recently, the lateral distribution of the measured radio 
signals has been investigated in more detail [4] . A function 
of the form 

e = e ° 6XP (~i^) ^ 

has been fitted to the measured field strength as a function 
of the distance to the shower axis Ro. The measured field 
strength as a function of the distance to the shower axis is 
presented in Fig. [I] for a typical event. The characteristic 
length Ro for this example is of order of 100 m. 

The investigations reveal that there are different types of 
air showers. Most of them have a lateral distribution which 
is characterized by a constant Ro ~ 100 to 200 m, like the 
example shown in Fig.[TJ On the other hand, there are few 
showers with relatively flat lateral distributions and cor- 
responding values for Ro as high as 1000 m or 1200 m. A 
closer look indicates that the steepness of the fall-off Ro 
depends on the mean distance to the shower axis R m ean 
and the zenith angle of the showers 0. The reconstructed 
scale parameter Ro is plotted as a function of the relation 
(1 - sin0)R 

mean i n Fig.O A correlation between the two 
quantities can be recognized. Large values for Ro are ob- 
tained for showers with large zenith angles and a small av- 
erage distance between the shower axis the antennas. But 
as not all events with small Rmean show a flattening — the 
reason is still unclear and further investigations with larger 
statistics are required. 

3. Polarization 

An important contribution towards the full understand- 
ing of the emission mechanism of the radio waves in air 
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Fig. 4. Ratio of the reconstructed cross-correlated beam values of the 
north-south and east-west polarization components as a function of 
azimuth angle (top) and the geomagnetic angle (bottom) |26| . 

showers is the investigation of the polarization character- 
istics of the measured radio signals [25]. The 30 LOPES 
antennas are set up now such that 15 register the east- west 
component and 15 antennas measure the north-south com- 
ponent of the electric field. 

The radio emission generated by the geo-synchrotron 
mechanism is expected to be highly linearly polarized. The 
signal is usually present in both polarization components 
(E-W and N-S). The signal strength depends on the geo- 
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Fig. 5. Example of the measured antenna signal as a function of time 
for a thunderstorm event 28 . 



magnetic angle and thus, on the shower azimuth and zenith 
angles [23]. The emission is expected to be polarized per- 
pendicular to the shower axis and the direction of the geo- 
magnetic field. The polarization characteristics can be de- 
scribed by a unit polarization vector v x B, where v is the 
direction of the shower axis and B the direction of the Earth 
magnetic field [26] . Recently, it has been suggested that in 
addition to the polarization characteristics also the abso- 
lute magnitude of the electric field is proportional to the 
Lorentz force v x B [27] . 

The arrival direction of the cosmic rays on the sky is 
shown for the registered events in Fig. 03 The color code 
represents the pulse heights of the cross-correlated beams, 
registered in east-west polarization (left) and north-south 
polarization (right). The boxes indicate the directions of 
the events with the strongest signals. The magnetic field 
in Karlsruhe has a zenith angle of 25° and an azimuth an- 
gle of 180°. It can be recognized that the showers with the 
strongest signals are registered perpendicular to the Earth 
magnetic field: from northern directions in east- west polar- 
ization direction and from west and east for the north-south 
polarized signal. 

To confirm the v x B behavior of the polarized signals 
the ratio of the amplitudes of the cross-correlated beams in 
east- west and north-south polarization is investigated [26] . 
The ratio of the north-south to east- west amplitudes is de- 
picted in Fig. 0] as a function of the azimuth angle (top) 
and the geomagnetic angle (bottom) . Mean values and the 
spread of the distributions are shown. A correlation be- 
tween the plotted quantities can be clearly recognized. The 
ratio of the north-south divided by the east-west compo- 
nent of the v x B vector of the measured showers exhibits 
the same behavior as a function of the azimuth angle and 
geomagnetic angle. This v x B behavior of the measured 
data is a strong indication for a geomagnetic origin of the 
radio emission in air showers. 
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Fig. 6. Sky map of the thunderstorm event shown in Fig. [5] [28 . 

4. Thunderstorm events 

The measured field strength of the radio emission of air 
showers depends on the (static) electric fields in the at- 
mosphere. The electric fields inside thunderstorm clouds, 
in particular within the convective region can reach peak 
values up to 100 kV/m. This leads to additional forces 
on the electrons and positrons in the air showers and, 
consequently, to an acceleration and deceleration of parts 
of the electromagnetic shower component (electrons and 
positrons, depending on the direction of the electric field). 
In turn, this yields to an amplification or reduction of the 
radio emission of air showers during thunderstorms. Such 
a behavior has been observed with LOPES [10J. To obtain 
reliable information about the detected air showers from 
the observations of radio emission requires to monitor the 
electric field in the atmosphere and to record the signa- 
tures of thunderstorms. Therefore, the signals form air 
showers during thunderstorms are studied in detail [28J. 

As an example, the antenna signal as a function of time 
measured during a thunderstorm is shown in Fig.[5j A sig- 
nal from an air shower is expected at the zero point of the 
time axis with an amplitude less than 0.1 V. Strong addi- 
tional signals caused by lightning can be recognized. The 
dashed lines mark a region with strong lightning signals, 
which is used to calculate a sky-map of the cross-correlated 
beam. The result is presented in Fig. [6] The map shows the 
whole sky with the zenith in the center and the horizon at 
the border. A strong signal is visible in south-western di- 
rection, marking the discharge region of the lightning. The 
signal extends to the horizon. This indicates that a cloud- 
ground discharge has been registered (and not a cloud-cloud 
discharge). The grating lobes of the antenna array extend 
over the whole sky. 

At present, studies are under way to detect a possible dis- 
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tortion of an event due to thunderstorm conditions directly 
from the measured data [28]. The idea is to find a possi- 
ble deviation of the polarization characteristics from the 
theoretically expected behavior. Such an indication would 
point towards a change in the emission process. 

5. New antenna types and self-trigger algorithms 

An important objective of LOPES is to develop and op- 
timize new antenna types to measure radio signals from 
extensive air showers in large installations, such as the 
Pierre Auger Observatory, see Sect.[U a sub-project named 
LOPES STAR . Most cosmic-ray radio detectors use (simple) 
dipole antennas, they are easy to assemble and are avail- 
able at relatively low costs. To avoid uncertainties related 
to dipoles new antenna types are developed, such as the log- 
arithmic periodic dipole antenna (LPDA) with two crossed 
polarization directions [30J. These antennas exhibit an al- 
most frequency independent directional diagram, antenna 
gain, and impedance. 

Another way to design wide-band directional antennas 
with dimensions much smaller than the LPDA is given by 
resistively loaded aperiodic antennas with internal losses 
[29] . They also have excellent wide-band properties as the 
resistor load dominates in comparison to the capacitive 
or inductive reactance. For a large-scale radio detector 
a crossed polarized short aperiodic loaded loop antenna 
(SALLA) has been developed. It has a diameter of only 
100 cm (120 cm in the latest version), a weight less than 
2 kg, and material costs of about 60 EUR. A photograph of 
the antenna and the circuit diagram are depicted in Fig.0 
The damping resistor guarantees a wide bandwidth. The 
E-plane directional diagram is shown in Fig. [7] (right). It 
features a wide main lobe towards zenith. The 3 dB beam 
width amounts to 150° and is about 50° wider than the 
directional pattern of the LPDA. 

Ultimate goal is to use these antennas in a self-triggered 
mode, i.e. the data acquisition is started directly from the 
detected radio signal [31J. A trigger algorithm has been 



developed and realized in FPGA-hardware. It provides RFI 
suppression by Fourier transforming the radio signal to the 
frequency domain in real time, eliminating mono frequent 
carriers and transforming it back to the time domain. Then 
a threshold is applied and events are selected which fulfill 
certain criteria to characterize the pulse shape. In a last 
step, coincidences between three neighboring antennas are 
required for a valid event. The real-time feasibility of this 
trigger mechanism was proven on prototype hardware. 

6. Outlook 

The investigations of radio signals from air showers with 
LOPES and CODALEMA are the basis for the application 
of these technique in large-scale experiments. The next step 
is to utilize the radio detection technique in large arrays, 
comprising several hundreds of antennas. 

LOPES serves as prototype for the detection of radio 
emission from air showers in LOFAR [32J. The Low Fre- 
quency Array (LOFAR) is a new digital radio observatory, 
presently under construction in the Netherlands and in Eu- 
rope [33]. It is designed as multi-sensor network to assist 
scientists in the fields of astronomy, geophysics, and agri- 
culture. Main focus of the astronomy community is to ob- 
serve the radio Universe in the frequency range from 30 
to 240 MHz. An objective of LOFAR is the detection of 
radio emission from particle cascades, originating from ex- 
tremely high-energy particles from outer space. Two main 
lines of research are followed: (i) the measurement of radio 
emission from extensive air showers, generated by inter- 
actions of high-energy cosmic rays in the atmosphere [34] 
and (ii) the detection of radio emission of particle cascades 
in the Moon, originating from ultra high-energy neutrinos 
and cosmic rays interacting with the lunar surface [35] . 

More than 40 stations with fields of relatively simple an- 
tennas will work together as digital radio interferometer. 
The antenna fields are distributed over several countries in 
Europe with a dense core in the Netherlands. The latter 
will have at least 18 stations in an area measuring roughly 
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2x3 km 2 . Each station will comprise 96 low band antennas, 
simple inverted V-shaped dipoles (like the LOPES anten- 
nas), operating in the frequency range from 30 to 80 MHz. 
Each antenna will have a dipole oriented in north-south 
and east-west directions, respectively. In addition, fields Fl 
of high-band antennas will cover the frequency range from 
110 to 240 MHz. For air shower observations the signals 
from the low band antennas are digitized and stored in a 
ring buffer (transient buffer board, TBB). For valid trig- 
gers the data are send to a central processing facility, based 
on an IBM Blue Gene supercomputer. First data from LO- 
FAR are expected in early 2010. 

Goal of LOFAR is to further push the development of ra- 
dio detection to be a new, independent way to measure the 
properties of air showers. Ultimate goal is to derive informa- 
tion about the primary, shower-inducing particle from the 
measurements, such as the particles energy, mass, direction 
and point of incidence. The curvature of the shower front 
has been investigated [36 37J. It could be shown that the 
radius of curvature measured at ground level from the radio 
observations is related to the distance to the shower maxi- 
mum. The depth of the shower maximum in the atmosphere 
Xmax is an important observable to determine the mass of 
the primary particle. The study indicates a resolution of 
X max from radio observations of order of 30 — 40 g/cm 2 . 

Another important application is the measurement of ra- 
dio emission from air showers at the Pierre Auger Obser- 
vatory with the Auger Engineering Radio Array (AERA) 
[38 39J. AERA is co- located with the infill array of the 
Pierre Auger Observatory and in the field of view of the 
high- altitude fluorescence telescopes (HEAT). This unique 
set-up will allow to register air showers simultaneously with 
three independent detection methods: radio waves, fluores- 
cence light, and particle detection in water Cerenkov detec- 
tors. AERA comprises about 150 antennas located on an 
area of about 20 km 2 and is designed to cover the energy 
range from 10 1T to 10 19 eV. It has three main science goals: 
(i) the thorough investigation of the radio emission from 
air showers and the analysis of the observed field strength 
as a function of various air shower parameters, such as en- 
ergy, angle of incidence, and distance to the shower axis, (ii) 
Exploration of the capability of the radio detection tech- 
nique and to study the feasibility as stand-alone technique 
to measure the properties of air showers, (iii) To measure 
the composition of cosmic rays in the energy range from 
10 17 to 10 19 eV. First data from AERA are expected before 
Summer 2010. 

The radio detection of air showers is a fast growing sub- 
discipline in astroparticle physics. With the new experi- 
ments starting operation exciting results are expected in 
the next few years. 
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1 The fields comprise 48 antennas in the Dutch stations and 96 in 
the European ones. 
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